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Abstract 

From the partial differential Calogero's (three-body) and Smorodinsky-Winternitz 
(superintegrable) Hamiltonians in two variables we separate the respective angular 
Schrodinger equations and study the possibilities of their "minimal" VT symmetric 
complexification. The simultaneous loss of the Hermiticity and solvability of the 
respective angular potentials V(<p) is compensated by their replacement by solvable, 
purely imaginary and piece- wise constant multiple wells Vo(<p). We demonstrate that 
the spectrum remains real and that it exhibits a rich "four series" structure in the 
double-well case. 
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1 Introduction 



One of the most important methods of solution of partial differential Schrodinger 
equations in D dimensions is an approximate [1] or exact [2] separation of variables 
in spherical coordinates. For example, once the potential is central, V(r) = V(\r\), 
the angular part of each wave function coincides with the well known D dimensional 
(hyper) spherical harmonics. 

In non-central case the explicit examples range from the confinement of a single 
(quasi)particle in a superintegrable Smorodinsky-Winternitz potential [3] up to the 
binding of A > 1 particles in the exactly solvable models of the Calogero-Moser- 
Sutherland type [4, 5]. Of course, the structure of the corresponding D dimensional 
generalizations of the (hyper) spherical harmonics becomes more complicated. Our 
present study was inspired precisely by this problem. 

In section 2 we shall introduce some basic notions and a few explicit more or less 
standard examples of the above-mentioned partial differential Schrodinger equations 
(cf. subsection 2.1 for more details). Then we shall specify our main purpose, 
namely, the analysis of the models in question in their appropriate non-Hermitian 
generalizations (a few remarks about motivation will be collected in subsection 2.2). 
In such a setting, subsection 2.2 will summarize the key technical aspects of our 
project. 

In section 3 we shall specify and explain in some detail the present implemen- 
tation of the Bender's and Boettcher's idea [6] of the weakening of the traditional 
Hermiticity to the mere VT symmetry. Subsection 3.1 will show that once we de- 
cided to study just the two-dimensional partial differential Schrodinger equations, 
their angular part degenerates to the mere single ordinary differential equation, 

^)=^ 2 ^), v?e(0,27r). (1) 

This simplifies many technical aspects of our problem which may be re-formulated 
as a confinement of a quasi-particle on a circle [7]. Unfortunately, in contrast to the 
standard models which use the Hermitian Hamiltonians and in which both the inter- 
pretation and the solution of the angular Schrodinger eq. (1) remain straightforward, 
we encounter a really serious difficulty here since the weakening of the Hermiticity 
makes our potentials V(ip) unsolvable. As we shall show below, they all acquire a 
multi-well shape with some 2M high barriers and deep minima all over the circle. 
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Recently we VT symmetrized the three-body Calogero model only via a very 
specific complexification of its angular part (1) which enabled us to find new solutions 
and to write them all down in closed form [8] . This admitted a preservation of many 
parallels between the Hermitian and non-Hermitian case. Still, one cannot assert that 
our recipe was fully satisfactory since, first of all, we did not provide any guarantee 
that the construction of the solutions in polynomial form does give us all the existing 
solutions. For this and other reasons we are going to suggest now an alternative VT 
symmetrization recipe. 

The details will be described in subsection 3.2 where we postulate that all our 
angular potentials possess just a suitable piece-wise constant shape which is, in addi- 
tion, purely imaginary. In this way, our choice will parallel many popular and purely 
imaginary VT symmetric models on the real line. Simultaneously, we merely lose 
a small part of the simplicity of the central models where the angular Schrodinger 
equation may be re-interpreted as a square well in quantum mechanics on a circle. 
For a more detailed explanation of such a point of view, cf. the first paragraph in 
subsection 2.1 below. 

By construction, our multiple imaginary square wells on a circle remain solvable 
by the matching method in principle. In practice, this is illustrated in section 4 
where we contemplate the simplest nontrivial two-dimensional problems of the above- 
mentioned partial differential form. In technical terms, this section deals with the 
angular Schrodinger equation defined on a circle and is a core of our present message. 
Its subsection 4.1 describes some details concerning the implementation of the current 
matching method, and the subsequent subsection 4.2 offers an ample illustration of 
this recipe and its results. In three Tables we sample the resulting spectra and discuss 
their apparent split in four sub-series. Some methodical aspects of the practical 
search for these energies are clarified with the assistance of several Figures. They 
document the reliability of the method and show that just keeping in mind a few 
simple rules concerning the (numerical) scaling, one may calculate the energies with 
any pre-determined precision. 

In our final section 5 we place all our results in a broader context, i.a., of the pop- 
ular double- well problems in subsection 5.1 and of their specific "angular" features 
in subsection 5.2. 
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2 Angular Schrodinger equations 



2.1 Hermitian non-central models 

Let us recollect the simplest possible two-dimensional situation where the angular 
Schrodinger equation (1) is just one-dimensional and defined on a circle (i.e., with 
the periodic boundary conditions). In all the central cases it contains just the trivial 
potential term, V((p) = for all (p G (0, 2ir). Hence, the well known (hyper)spherical 

solution of eq. (1) is immediate, i/) m ((p) = exp imip, £ = m — 0, 1, 

In all the less trivial but still fully separable and exactly solvable models with non- 
central potentials, the angular Schrodinger equation (1) remains exactly solvable by 
assumption. One may recollect one of the oldest and simplest illustrative examples 
of this type, viz., the Smorodinnsky-Winternitz "superintegrable" oscillator in two 
dimensions [3], 



d 2 d 2 v2 v2 G G 
+ X 2 + Y 2 + — + 



$(X,Y)=S${X,Y). (2) 



OX 2 dY 2 X 2 Y 2 

After a routine elimination of the uninteresting motion of the centre of mass, also the 
partial differential representation of the Calogero's popular three-body bound-state 
model [4] has the very similar form 

d 2 d 2 



dX 2 dY 2 



X 2 + Y 2 + 



$(X,Y)=£$(X,Y). (3) 

2X 2 2{X-V3Y) 2 2{X + V3Y) 2 _ 
Due to the presence of the repulsive barriers at G > and g > 0, the angular parts 
of both these equations cease to be defined on the corresponding hypersphere (i.e., 
on the circle). Thus, once we distinguish the above two models by their respective 
"label numbers" M = M Winternitz = 2 and M = M Ca iogero = 3, we may summarize 
that 

(a) the domain of our innovated harmonics ip{(p) shrinks to a segment tp G (0, tt/M) 
of the full circle, 

(b) at both the ends of this interval, the Dirichlet boundary conditions must be 
imposed, due to the presence of a strong singularity there, 

(c) both the models (with the elementary Poschl- Teller angular potentials) still 
remain exactly solvable, as expected, in terms of Jacobi polynomials [9]. 
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All the three comments (a) - (c) remain valid also for an oversimplified "toy" model 



<9 2 <9 2 G 
° _ j± I y 2 i y 2 i y 

dx 2 dY 2 x 2 



<S>(X,Y) =S${X,Y). 



(4) 



Its "index" M = M toy = 1 is "minimal" and the model is popular as a not entirely 
trivial methodical guide [4, 8]. It will play the same guiding role in what follows. 



2.2 A motivation for transition to non-Hermit ian Hamilto- 
nians 

In the recent literature, an increasing attention is being paid to non-Hermitian Hamil- 
tonians with real spectra. Mostly they are characterized by their parity-times-time- 
reversal symmetry, H = VTHVT [6, 9, 10] or, in a more formal setting, by their 
specific parity-pseudo-Hermiticity property H = VH^V [11] [13]. Via analytic 
continuation this tendency seems to offer a nontrivial insight in the solutions [14]. 

An important phenomenological appeal of the similar studies is related to the 
observation that the complex conjugation operator T mimics the time reversal [15, 
16]. The early stages of development of this subject were co-motivated by its natural 
emergence within perturbation studies of various models in quantum mechanics and 
field theory. People revealed, for example, that even the peculiar and "strongly 
non-Hermitian" imaginary cubic anharmonic oscillator may have a real and discrete 
spectrum bounded from below [17, 18]. A few years later, an introduction of an 
"unstable" quartic anharmonic oscillator proved mathematically relevant in having 
facilitated the perturbative description of some apparently non-perturbative double- 
well models (e.g., of H^ DW \g) = p 2 + x 2 (l — g x 2 ) 2 + cx in ref. [12]). The weakening 
of the Hermiticity to mere pseudo-Hermiticity was even found to make the non- 
Hermitian quartic oscillators exactly solvable at certain couplings and energies [19]. 

The quick progress in the field has led to the sophisticated perturbative anal- 
yses of some non-Hermitian VT symmetric models in more dimensions [20]. The 
"analytic continuation" activities intensify since the Hermitian non-central partial 
differential Schrodinger equations are attractive by the variability of their possible 
physical interpretations. They might prove relevant in the models which range from 
a single particle moving in a non-central D dimensional potential [21] up to a full- 
fledged A body motion and spectra [22] . Still, the real impact of their non-Hermitian 
analytic continuations is not yet fully predictable at present. 
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In contrast to the majority of contemporary publications which study the VT 
symmetric models in one dimension, one also has no clear guide to the purely tech- 
nical aspects of the VT symmetrization of the models in more dimensions. We 
feel motivated by this mathematical challenge and intend to pay a deeper attention 
to the specific subset of the VT symmetrization techniques which is based on an 
"incomplete" complexification involving just the "angular" part of the problem. 

Without any detour to perturbation methods we shall study some specific prop- 
erties of a class of the "angular" Schrodinger equations in a way encouraged by 
our experience with the comparatively easy tractability and insight provided by the 
piece-wise constant potentials [23]. 

2.3 Complexifications of the angular Schrodinger equations 

Within the framework of the so called VT symmetric quantum mechanics [16] one 
may complexify the model (3) of Calogero [8, 24] as well as its spiked and superinte- 
grable predecessors (2) [25] by the complexification of its single angular coordinate. 
Such an approach remains sufficiently transparent. One of its very friendly features 
is that it leads to a perceivable strengthening of the analogies between many above- 
mentioned non-central models [with g ^ in eq. (3) or with G ^ in eq. (2)] and 
their central, regular special cases (with g = or G = 0, respectively). 

In what follows, we just intend to study these analogies in more detail. In a more 
explicit formulation of this purpose, our analysis will start form the weakening of 
the Hermiticity to the mere VT symmetry of the Hamiltonian. This opens a way 
towards a suitable complex deformation of the integration path [6], especially in the 
vicinity of the singular points ifk = nk/M, k = 1,2, ...,2M. In this sense one 
should 

(A) re-extend the domain of the complexified harmonics ip(<p) to the full circle, 

(B) return to the standard periodic boundary conditions in order to keep our wave 
functions unambiguously defined, 

(C) re-analyze the solutions of eq. (1) in complex domain. 

Unfortunately, none of the existing papers did solve all the three problems (A) - (C) 
at once. In particular, paper [26] studied the solutions of eq. (1) on a single, isolated 
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complex segment, while a return to the full circle has been facilitated in refs. [8, 25] 
by the use of the "non-periodic" , Dirichlet boundary conditions. 

An incompleteness of the latter results formed a key inspiration of our present 
paper. In essence, it will circumvent some of the main technical obstacles via a 
simple-minded replacement of the complexified Poschl- Teller potentials in eq. (1) by 
their square-well-like solvable approximants. 

3 Toy models on a circle 

3.1 A constant complex shift of the angle <p 

Let us consider one of the above-listed angular Schrodinger equations (1) with the 
Poschl- Teller M— dependent potential. One of its simplest complexifications on a 
circular domain will be then determined by a mere constant complex shift (plus a 
comfortable re-scaling) of the angular variable, 



Whenever the imaginary shift is very large, a»l, the Poschl- Teller potential itself 
becomes very smooth and "almost forgets" about all its singularities along the real 
line, 



This potential is non-Hermitian and manifestly VT symmetric (with the current 
definition of the complex conjugation T and parity Vs = —s). In generic case, its 
"energies" [or rather eigenvalues £ 2 in eq. (1)] may be expected real [16]. In addition, 
its wave functions may be comfortably expressed in terms of Bessel functions [27]. 

Unfortunately, the above solutions only become exact in the limit a — > oo where 
the potential (6) itself vanishes. Obviously, all the information about the original 
angular potential is lost and nothing new is gained in the limit a — > oo since, formally, 
one just returns to hyperspherical harmonics. Vice versa, trying to preserve the 
specific, 2M— well character of our family of the Poschl- Teller angular potentials 
y( PT ) we must return to the domain of the finite shifts a. In what follows, we intend 
to make just the first steps in this direction. 



<P = V( s ) = -ia + -s, 



a>0, se(-2,2). 



(5) 



V( PT \ip) ~ l/sm 2 (Mip) = e~ 2Ma ■ - 



cos ttMs + i sin irMs + O (e' 2Ma ) ] . (6) 
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3.2 Piece-wise constant angular potentials 

Basically, we shall preserve the multi-well character of our potentials, trying merely 
to deform the local details of their global shape. In a more detailed specification 
of our "toy-model" we were guided by refs. [28] and [29]. In the former one, the 
complex and VT— symmetric potentials with elementary trigonometric shapes re- 
sembling eq. (6) have been analyzed in periodic case. In accord with expectations, 
the corresponding band spectrum proved real but in contrast to the conventional 
real periodic potentials, it appeared to exhibit certain significant qualitatively new 
features. Although, due to the unsolvable character of the potentials, all these obser- 
vations were originally based just on the careful numerical analysis and higher-order 
WKB techniques, their confirmation has been delivered in the subsequent short note 
[29]. There, the presumably not too relevant real part of the "realistic" potential 
V has been replaced by a constant (zero). Moreover, the exact solvability of the 
new schematic model has been achieved by an extremely drastic reduction of all the 
remaining (i.e., purely imaginary) part of the potential V to the mere sequence of 
delta functions (known as the Kronig- Penney model) with purely imaginary (i.e., still 
VT— symmetric) strengths mimicking the presence of the infinitely many minima in 
the original V. 

Here we shall proceed in a very similar manner. Firstly, we shall ignore again 
the real part of the potential as not too relevant. Besides all the above-mentioned 
indirect motivation of such a technical simplification, we also directly verified its 
consistency in ref. [30] where a "minimal" non-periodic double-well version of the 
Kroning-Penney potential has been used. 

Secondly, we shall keep in mind that the most important feature of the remaining, 
purely imaginary part of our potentials (6) lies in the presence of their precisely 2M 
separate extremes. Thus, a transition to the solvable model will be achieved here via 
their piece-wise constant approximation. In this manner, once we denote h = 1/M, 
we arrive at the angular Schrodinger eq. (1) in the more or less unique form 



ds 2 



+ V {M \s) ^ M \s)=E^ M \s), 



se(-2,2) 



(7) 



with the very specific class of the purely imaginary toy potentials 




(8) 
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where 



A< + ' = (-2, -2 + h) U(-2 + 2/i, -2 + 3/i) U • • ■ U( 2 ~2h,2-h), 

and 

A (_) = (-2 + /i, -2 + 2/i) U . . . U( 2 - 3/i, 2 - 2/i) |J(2 - /i, 2) . 

For illustration, Figure 1 shows the shape of such a potential defined over the circle 
at the Calogero's choice of M — 3. 

In what follows we shall consider our potentials (8) and concentrate on their first 
nontrivial special case with M — 1, 

1 -«Z se (-1,0)11(1,2) 

This means that a combination of the VT symmetry requirement with the double- 
well character of the potential makes our M — 1 model defined on a circle of the 
convenient radius R = 2/n. 



4 Solutions 

4.1 The sample of the matching conditions at M = 1 

The explicit solutions of our toy example (9) will be constructed by the matching 
method. For this purpose, it makes sense to denote the left and right intervals (—2, 0) 
and (0, 2) by the respective subscripts l and r. The subintervals which lie near or 
far from the origin at s = will be assigned another subscript jv or F , respectively. 
In this notation we have the natural ansatz 

ipjix) = Aj sin KjX + Bj cos KjX, j = FL, NL, NR, FR (10) 

with 

Kpi, = knr = h = s — it, knl = kfr = k* = s + it, s > 0, t > 0. (11) 

This should be interpreted as a mere suitable reparametrization of the desired spec- 
trum of "energies" 

E = f = s 2 -t 2 , 2st = Z>0. (12) 
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Now its is easy to require the matching of the logarithmic derivatives at the discon- 
tinuity at x = — 1, 



^ FL (-1) = Vatl(-I), ^k(-l) = ^l(-I) , (13) 

and, mutatis mutandis, at the other two "internal" boundaries at x = and x — 1. 
The points rr = 2 and rr = — 2 must be matched together in a way which reflects the 
circular symmetry of the problem, 

if> FR (2) = VW-2), if/ FR (2) = if/ FL (-2) . (14) 

We may conclude that the elementary insertions lead to the linear system of eight 
homogeneous equations for all the eight coefficients Aj and Bj (which may be ar- 
bitrarily normalized) . This means that the secular determinant of the system must 
vanish, 

det Q(s, t) = . (15) 

After one inserts, say, s = s(t) = Z/2t from eq. (12), the insertion of all the non- 
vanishing matrix elements leads to the closed and compact form F(t) = of eq. (15) 
with 

Qu = ~Q*i3 = -Q55 = Qh = - sin k (16) 

Qu = -Qu = Qk = -Q™ = = -Ql 8 = cos K ? ( 17 ) 

Q21 = -Q23 = Q*43 = = Qd5 = -Q*67 = KCOSK (18) 

Q22 = ~Q*24 = ~Qm = Q*68 = KSinK, (19) 

Q71 = Q*77 = -sin2/t, Q 72 = -Q* 78 = cos 2k, (20) 

Qsi = ~Q*87 = K cos2/t, Q 82 = Q* 88 = /tsin2/t (21) 

and 

F(t) = -8 1 4 sinh 6 t sin 6 Q^j +...-— cosh 6 t cos 6 (J^j . (22) 

We did not display it here in full detail because with all its 30 separate real trigono- 
metric terms, this explicit formula is quite lengthy and its detailed inspection seems 
hardly instructive at all. For this reason, we are rather going to display here a few 
characteristic graphs of the function F(t) instead. 
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4.2 The sample of a graphical determination of the spectrum 

One should start in a weak coupling regime. Choosing the very small Z — 1/10 we 
get the sequence t = 0.2219819562, ti = 0.03467067057, . . . of the real roots of F(t). 
Due to the fact that the function F(t) itself has the elementary form, the numerical 
values of these roots may easily and quickly be evaluated with an arbitrary precision. 
These values then determine the low-lying energies E = 0.00153255, E\ = 2.078577 
etc. 

An increase of the coupling to Z = 1 shifts the lowest roots (to = 0.6564195696 
etc) and energies (E = 0.1493123386 etc) upwards. No changes in the pattern are 
observed, and the more detailed description of the Z = 1 results is given here in the 
three Tables 1-3. 

These sample results reveal and demonstrate a certain quadruple periodicity of 
the energy levels. This merely reflects the simplicity and regularity of the secular 
determinant F(t) itself. An inspection of an overall shape of this function clarifies 
the structure of the spectrum in a way which complements very well the explicit 
Tables. Thus, we may see in Figure 2 that beyond its rightmost or "maximal" root 
t ~ 0.66, the curve F(t) remains negative and decreases towards — oo very quickly. 

In a way, the ground-state root to has a certain "exceptional" status in having 
no natural "neighbor". For all the rest of the spectrum we observe that the roots 
may be organized as certain doublets (t 2 k-i, t 2 k) numbered by k — 1, 2, . . .. With an 
increase of the index k, these partners in the doublet lie closer and closer to each 
other. This is illustrated in Figure 3 where an onset of the existence of the doublets 
is clearly visible near t\ = 0.3934710177 and where we also notice a rapid growth 
of the maxima of F(t) which makes the numerical search for its zeros slightly more 
complicated. In particular, it is necessary to keep the necessary re-scalings of F(t) 
during its graphical analysis under control. 

For this purpose we used an algorithm based on an interplay between the plot- 
ting and root searching steps. For example, the doubts over Figure 3 concerning 
the existence of the root near t ~ 0.16 are easily settled by the magnification of 
this subinterval in Figure 4. In this way we were able to detect and scale-out the 
comparatively quick increase of the density of the roots t n > with the increase of 
their subscripts. 

A new phenomenon emerged in the domain of n ^> 1, which may be even noticed 
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in the Figure 4 itself. One detects a systematically growing difference between the be- 
haviour of the doublets (^fc-i, ^2k) at their even and odd indices k. The phenomenon 
is even better visible after we move to the energies themselves. Thus, a glimpse in 
Table 1 verifies immediately that only at the even k = 2m the quasi-degeneracy of 
the pairs of roots (t 2 k~i,t2k) has a real impact. Only in this case it implies also the 
similar quasi-degeneracy effect for the related energies, E im _i E im , m >> 1. 

The comparatively bigger distance between £2*1-1 and tik at the odd k = 2m + 1 
makes the numerical determination of these pairs less difficult. Vice versa, the quick 
decrease of distance between these roots at the even k = 2m may be perceived, 
especially in combination with the growth of the maxima of F(t) at the smaller 
t — > 0, as the key obstacle encountered during the precise numerical localization of 
the energies. A useful tool of a disentanglement of these two obstacles has been 
found in the introduction of an auxiliary function L(t) = log[|F(t)|], the practical 
numerical picture of which is sampled in Figure 5. 

Two of its properties have to be emphasized. Firstly, the absolute-value-part of 
the definition of L(t) guarantees that the logarithm remains well defined and that 
the changes of sign of the original functions F(t) on a very small interval become, 
empirically, "magnified" . The emergence of every quasi-degenerate pair of zeros will 
be reflected by a very well visible "artificial bump" in L(t) in place of a hardly visible 
doublet of nodes in an improperly rescaled graph of the unmodified function F(t) 
(compare, e.g., Figures 4 and 5). 

Secondly, the sequence of the maxima of L(t) will not grow that quickly with 
the decrease of t <C 1. Hence, the growing density of the roots is better represented 
by the logarithm, in spite of its, in principle, unlimited decrease near the zeros of 
F(t). Empirically, the latter effect proves negligible in numerical practice, in the way 
illustrated very persuasively by Figure 5. 

We may summarize that the main experience gained during the graphical solution 
of our sample secular eq. (15) lies in the necessity of an appropriate selection of the 
interval in which the separate roots t n are being sought. Of course, the "good" size 
of these intervals will depend on the coupling Z so that in general, one should never 
forget to consult the graphs of F(t) and L(i) . 
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5 Discussion 



5.1 Multiple- well problems in non-Hermit ian framework 

We may summarize that a family of certain purely imaginary potentials was con- 
sidered on a circle. Their piece-wise constant shape was chosen to mimic the multi- 
well angular part of certain not yet fully understood pseudo-Hermitian versions of 
Calogero-Moser and Smorodinsky-Winternitz Hamiltonians with real energies. Vice 
versa, the reality of our spectrum reflects a "minimal" square-well-type VI sym- 
metrization of the initial many-body or more-dimensional separable Hamiltonian. A 
thorough numerical analysis of the spectrum was performed for the first nontrivial 
special case with double well potential. It has been shown to exhibit a rich "four 
series" structure. 

These results fit well a broader perspective. Forty years ago, perturbation analysis 
of the sing le-well oscillator H (sw \g) = p 2 + x 2 + g 2 x 4 revealed that the most natural 
interpretation of its "physical" spectrum {e n (g)} emerges via a complexification of 
its coupling g [31]. All the apparently uncorrelated numbers e n (g) were shown to 
coincide with the values of a single analytic function e(g), evaluated on the different 
Riemannian sheets. In this interpretation, the main quantum number n = 0, 1, . . . 
merely selects the different Riemannian sheets in a certain specific, very transparent 
and intuitively appealing manner. One may conclude that the complexification of 
the parameters is one of the best ways towards an understanding of the spectra, 
especially when they exhibit a strongly anharmonic double well shape [32] since, in 
these cases, there may exist a close connection between the single and double wells 
mediated by Fourier transformation and analytic continuation [12]. 

A priori one might expect that the high degree of an internal symmetry of the 
superintegrable equations (which implies their separability and exact solvability in 
several systems of coordinates [21]) might mediate the facilitated technical feasibil- 
ity and, hence, preferrence of their systematic VT symmetrization. For this rea- 
son, our present recipe might be re-read as a direct "brute-force" introduction of a 
Hermiticity-breaking term. Such a modification of the superintegrable Schrodinger 
equations by an additional potential might offer a new insight in the problem of the 
reality of the spectra and, in particular, in the (sub)problem of the mechanisms of 
its possible (and, in principle, admissible and often quite easy) breakdown. 



12 



5.2 Summary and outlook 

The relationship of the symmetry of the Hamiltonian to the reality of its spec- 
trum was illustrated here via a few most elementary non-Hermitian solvable models 
in two dimensions. In the similar examples, the choice of the suitable VT symmetric 
anharmonic terms is usually dictated by the purely technical requirements of sim- 
plicity and solvability. In this sense, one encounters serious difficulties when moving 
to two and more dimensions since only the one dimensional Schrodinger equations 
abound with the solvable potentials. To a successful VT symmetric generalization of 
these models to more dimensions, the key contribution of our present short note may 
be seen in its proposal and verification of a combination of the matching and analytic 
techniques. We showed that in spite of the violation of the superintegrability caused 
by our introduction of the step-like imaginary potentials, we were still able to study 
our models by the essentially non-numerical technique. 

The angular potentials we considered were drastically simplified to a piece-wise 
constant purely imaginary field which manifestly destroys the Hermiticity of the 
Hamiltonian. Our extreme simplification of the generic situation offered the picture 
where the barriers become penetrable. This changes the boundary conditions to 
periodic ones [14]. The related analysis proved extremely illuminating and revealed 
the existence of the new and hardly expected types of a regularity in the spectrum. 
We may conclude that our choice of the schematic toy model (9) was fortunate, 
indeed. 

• Qualitatively, the energies in the spectrum seem to be classified as lying in the 
four separate series. 

• The intrinsic double-well character of the model seems to be responsible for 
the emergence of the almost degenerate energy doublets. 

• One encounters the quasi-degeneracy of the energies not only in the low energy 
domain but, quite unexpectedly, all over the energy scale. 

There exist several open questions which remained unanswered. Firstly, a deeper 
understanding of many empirical observations is missing. The quadruple-series spec- 
trum in Table 1 as re-emphasized in Tables 2 and 3 might motivate a continuing 
analysis of the more-dimensional systems. 
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Secondly, there is no doubt that the structure of the spectrum in our present 
VT symmetric model is richer than in its Hermitian predecessors. Although all 
the intuitive and/or elementary interpretations of the "angular-type" motion over 
the circle remain similar, one feels tempted to emphasize the differences. In this 
sense, one notices, first of all, that our present model re-confirms the idea and/or 
existence of certain "quantum beats" as conjectured on the basis of a different model 
in ref. [30]. 

Thirdly, one has to re-emphasize that in contrast to the quasi-degeneracy as 
observed in the Hermitian multiple wells, the present, PT-symmetric analogue of 
this phenomenon is "periodic" and transcedes the standard limitation of the quasi- 
degeneracy to the mere low-lying spectrum. 

Finally, we might re-open the discussion which has been initiated in ref. [14] . An 
ambiguity of a return to the Hermitian case has been revealed there, and interpreted 
as a transformation which gives rise to the unexpected emergence of some entirely 
new models. Of course, once we started here from a separation of a circle in some 2M 
quasi-independent sectors, similar outcome might be expected in principle. Indeed, 
our example just deals with the parts (quadrants) of a bigger phase space. As 
long as this space was originally decomposed in a few invariant subspaces, its non- 
Hermitization might offer a recipe how one could open a tunnel through the originally 
fully impenetrable barriers. 
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Table captions 



Table 1. The low-lying energy spectrum at Z = 1. 

Table 2. Two natural subsets of the even energy levels of Table 1. 

Table 3. Two natural subsets of the odd energy levels of Table 1. 



Figure captions 

Figure 1. Imaginary part of a Calogerian six- well angular potential in its rectan- 
gular approximation. 

Figure 2. Monotonous decrease of the secular determinant F(t) beyond its 
ground-state root t 0.656 at Z — 1. 

Figure 3. Graphical search for the roots of F(t). 

Figure 4. A magnification of Figure 3 near t = 0.16. 

Figure 5. Logarithmic re-scaling L(t) = log[|F(t)|] and an improved graphical 
trapping of quasi- degeneracy. 
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Table 1: The low- lying energy spectrum at Z = 1 



level 


parameter 


energy 


difference A n 


n 


tn 


E n 


E n — E n _i 





0.656 


0.149312 


- 


1 


0.393 


1.459965 


1.319 


2 


0.266 


3.464686 


2.004 


3 


0.1596 


9.792771 


6.328 


4 


0.1587 


9.895111 


0.102 


5 


0.111 


20.127356 


10.232 


6 


0.101 


24.273237 


4.146 


7 


0.07959 


39.459389 


15.186 


8 


0.07956 


39.484770 


0.025 


9 


0.0651 


58.954565 


19.469 


10 


0.0623 


64.410749 


5.456 


11 


0.053053 


88.817991 


24.397 


12 


0.053050 


88.829258 


0.011 


13 


0.0461 


117.64814 


28.819 


14 


0.0449 


124.15474 


6.507 


15 


0.039789 


157.90892 


33.754 


16 


0.039788 


157.91526 


0.006 


17 


0.0357 


196.15402 


38.239 



Table 2: Two natural subsets of the even energy levels of Table 1 



level 


second difference 


level 


second difference 


n 


A n A n _4 


n 


A ra — A ra _4 


6 


2.142 


8 


-0.77 T 


10 


1.310 


12 


-0.014 


14 


1.051 


16 


-0.0049 



Table 3: Two natural subsets of the odd energy levels of Table 1 



level 


second difference 


third difference 


n 


A n - A n _ 2 


A n - 2A n _ 4 + A n _ 8 


3 


5.018 




5 


3.904 


-1.114 


7 


4.954 


+ 1.050 


9 


4.283 


-0.671 


11 


4.938 


+0.655 


13 


4.422 


-0.516 


15 


4.935 


+0.513 


17 


4.485 


-0.450 
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